Quantum Hall effect in a one-dimensional dynamical system 
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We construct a periodically time-dependent Hamiltonian with a phase transition in the quantum 
Hall universality class. One spatial dimension can be eliminated by introducing a second incom- 
mensurate driving frequency, so that we can study the quantum Hall effect in a one-dimensional 
(ID) system. This reduction to ID is very efficient computationally and would make it possible to 
perform experiments on the 2D quantum Hall effect using cold atoms in a ID optical lattice. 
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I. INTRODUCTION 

The disorder-induced localization-delocalization tran- 
sition in the quantum Hall effect is the oldest and best- 
known example of a topological phase transition*^ The 
transition is called topological because it is associated 
with a change in a topological invariant, the Chern num- 
ber, which counts the number of edge states and the 
quantized value of the Hall conductance^ Since there 
is still no analytical theory for the quantum Hall phase 
transition, computer simulations are needed to calculate 
the scaling law and critical exponent associated with the 
diverging localization length at the transition. The two- 
dimensional (2D) network model of Chalker and Cod- 
dington has been the primary tool for these studies for 
more than two decades 

In this paper we introduce an alternative one- 
dimensional (ID) model of the quantum Hall phase tran- 
sition. The model is stroboscopic, with a Hamiltonian 
that is driven quasipcriodically with two incommensu- 
rate driving frequencies. It is a variation on the quantum 
kicked rotator^— used to study the 3D Anderson metal- 
insulator transition of atomic matter waves in a ID opti- 
cal lattice*^— Stroboscopic models of quantum phase 
transitions have received much attention recently^ - — 
but the dimensional reduction considered here was not 
yet explored. 

Usually the quantum Hall effect is due to the quantiza- 
tion of cyclotron orbits in Landau levels. It is possible to 
simulate a Lorentz force acting on neutral atoms in a 2D 
optical lattice^ - — but in a ID lattice we need a quantum 
Hall effect without Landau levels^ This socalled quan- 
tum anomalous Hall effect appears in the Qi-Wu-Zhang 
(QWZ) model of a spin-1/2 coupled to orbit and to a 
uniform magnetization. While the topological invariant 
in this model takes on only the three values 0,±1, the 
phase transitions are in the same Z universality class as 
the usual quantum Hall effect. 

In the next two sections we formulate the stroboscopic 
model of the quantum Hall effect, first in 2D (Sec.|H| and 
then reduced to ID (Sec. IIIip . We obtain the model by 
starting from the QWZ Hamiltonian, but we also show 
how it is related to the quantum kicked rotator (upon 



exchange of position and momentum). 

In Sec. IIVI we perform numerical simulations of the 
spreading of a ID wave packet to identify the localization- 
delocalization transitions. While the translationally in- 
variant QWZ model has three quantum Hall transitions, 
we find four transitions because one is split by disorder. 
We verify one-parameter scaling of the time-dependent 
diffusion coefficient and calculate the critical exponent. 
The result is consistent with the most accurate value ob- 
tained from the Chalker-Coddington models 

To further support that these are topological phase 
transitions in the quantum Hall universality class we cal- 
culate the Hall conductance as well as the Z topological 
invariant in Sec. El We conclude by discussing the pos- 
sibilities for the realization of the quantum Hall effect in 
a ID optical lattice. 



II. FORMULATION OF THE 2D 
STROBOSCOPIC MODEL 

A. Quantum anomalous Hall effect 

In this subsection we summarize the QWZ model^ of 
the quantum anomalous Hall effect, on which we base the 
stroboscopic model described in the next subsection. 

The QWZ model describes two spin bands of a mag- 
netic insulator on a two-dimensional (2D) square lattice. 
The crystal momentum p = (p\,P2) varies over the Bril- 
louin zone —Trh/a < Pi,P2 < nh/a. The Hamiltonian has 
the form u ■ a, with a = (cr x ,a y ,a z ) a vector of Pauli 
matrices and 

u(p) = K(smp 1 ,smp 2 ,/3\p - cospi - cosp 2 ]) (2.1) 

a momentum-dependent vector that couples the spin 
bands. (We have set H and a both equal to unity.) The 
dispersion relation is E±(p) = ±u(p), with u = \u\ the 
norm of the vector u. We fix the Fermi level at zero, in 
the middle of the energy gap. 

Eq. (|2.ip contains three parameters, K, (3, [i. The pa- 
rameter K sets the strength of the spin-orbit coupling. 
Time-reversal symmetry is broken by a nonzero /3, repre- 
senting a magnetization perpendicular to the 2D plane. 
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The Hall conductance Gh is quantized at ?3 
r sign/3^ if \fi\ < 2, 



G H = T x 







if |/z| > 2. 



(2.2) 



This quantum Hall effect is called "anomalous" , because 
it does not originate from Landau level quantization. 

The value of Gh is a topological invariant^ meaning 
that it is insensitive to variations of the Hamiltonian that 
do not close the energy gap. Since the gap can only close 
if u(p) vanishes for some p, the Hamiltonian 



Ho(p) = T{u) u ■ er 



(2.3) 



has the same quantized Hall conductance (|2.2[) if the 
function T{u) is positive definite. We will make use of 
this freedom in order to flatten the spin bands, by choos- 
ing a function T(u) which decays for large u. 

The Hamiltonian Ho(p) describes a clean system. The 
effects of electrostatic disorder are included by adding the 
scalar potential V(x). The 2D coordinate x = {x\,X2) is 
measured in units of a, while momentum p = (pi,p2) 
is measured in units of H/a, so their commutator is 



B. Stroboscopic Hamiltonian 

This completes the description of the QWZ model. We 
now introduce a periodic time dependence by multiplying 
Hq with the stroboscopic function r 8(t — tit), while 
keeping the scalar potential time-independent. We thus 
arrive at the stroboscopic Hamiltonian 

oo 

U(t) = V(x)+H (p) S(t-n), (2.4) 

n— — oo 

where we have set the period r equal to unity. 

For the choice of T{u) and V(x) we are guided by the 
tight-binding representation given in App. [X] We use 



T(u) 



2 arctan u 



(2.5) 



which has a tight-binding representation with nearest- 
neighbor hopping. For the scalar potential V(x) we take 
a separable form, 



(2.6) 



with Viixi) a low-order polynomial in Xi. Such a simple 
potential produces quasi-random on-site disorder in the 
tight-binding representation. 



C. Relation to quantum kicked rotator 

The quantum kicked rotator is a particle moving freely 
along a circle, with moment of inertia /, being subjected 




FIG. 1: Momentum dependence of the quasi-energy (|2.9|l for 
zero disorder potential, calculated from Eqs. (|2,l|l and (|2.5p 
for K = 2, f3 = 0.8, /i = 1.9. At the center of the Brillouin 
zone the Dirac cone emerges, which will be fully formed when 
the gap closes at /j, — 2. 



periodically (with period r) to a kick whose strength de- 
pends oc cos 9 on the angular coordinate 9. The quantum 
mechanical Hamiltonian isL2i 



H(t) 



2lW 



KI 



cos 9 s (t 



(2.7) 



n— — oo 



The stroboscopic Hamiltonian (|2.4[) has the same general 
form, upon substitution of 6 i— > p, with the extension 
from ID to 2D and with the addition of a spin degree of 
freedom in the kicking term. 

A ID spinfull kicked rotator has been used to 
study the effects of spin-orbit coupling on quantum 
localizationi^r— (Because in the kicked rotator the vari- 
able which localizes is momentum rather than position, 
one speaks of dynamical localization.) In these ID stud- 
ies there was only a topologically trivial phase, while — 
as we shall see — the present 2D model exhibits a topo- 
logical phase transition. 



D. Floquet operator 

The evolution ^(t + 1) = J-^(t) of the wave function 
^f(t) over one period is described by the Floquet operator 
T. Integration of the Schrodinger equation id^ /dt = 
T-L(t)^(t) gives the Floquet operator as the product 



_ e -iH {p) e -iV(id p ) 



T = e 



(2.8) 



with id p = id /dp the position operator x in momentum 
representation. 

The eigenvalues of the unitary operator J- arc phase 
factors e~ l£ . The phase shift e € [— ir, ir) plays the role of 
energy (in units of fo/r), and is therefore called a quasi- 
energy. For V = the quasi-energy is an eigenvalue of 
H , hence 



e = ±uT{u), for V = 0. 



(2.9) 
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The p-dependence of the two bands is plotted in Fig. [T] 
The emerging Dirac cone is clearly visible. Away from 
the cone the bands are quite flat, which is a convenient 
feature of our choice (|2.5p of T(u). 

More generally, for nonzero V, the 27r-periodicity of 
Ho(p) implies that the eigenstates 



**(p) = e- lpq Xq {p) 



(2.10) 



of T are labeled by a Bloch vector q in the Brillouin zone 
— 7r < <7i,<72 < 7T- The function Xq(p) is a 27r-periodic 
cigenstate of 

T q = e - lHa[p) e-' lV{ld f +q) . (2.11) 

A convenient basis for the functions Xq(p) is formed by 
the eigenfunctions exp(— im-p) of x. The 27r-periodicity 
of Xq requires that the vector m = (7711,7712) contains 
integers. 



III. MAPPING ONTO A ID MODEL 

The quantum kicked rotator in d dimensions can be 
simulated in one single dimension by means of d incom- 
mensurate driving frequencies We apply this dimen- 
sional reduction to our stroboscopic model of the quan- 
tum Hall effect. 

We take a linear potential in the variable X2, 



V(x) = Vi(xi) - u>x 2 , 



(3.1) 



with uj/2tt an irrational number s (0, 1). During one pe- 
riod the momentum P2 is incremented to P2 + w (modulo 
27r), so u is an incommensurate driving frequency. An 
initial state 

tf(pi,Pa,* = 0) = i/)(pi,t = 0)S( P2 - a). (3.2) 
evolves as 

*(pi,pa,t) = e- tHn< - Pl ' ut+ ^e- lV ^ Xl ^( Pl ,t - 1) 

x 5{p 2 -wt-a). (3.3) 

We may therefore replace the 2D dynamics by a ID dy- 
namics with a time-dependent Floquct operator 

F{t) = e -iffo(pi,u<t+a) e -iV 1 (ia J , 1 ) ) ( 34 ) 
T q (t) = e -i#o(pi,«t+a) e -iVi(ia J , 1 +?)_ (35) 

This reduction from two dimensions to one dimension 
greatly simplifies the numerical simulation of the quan- 
tum Hall effect. 

For the potential in the remaining dimension we take 
a quadratic form, 

V 1 {x 1 ) = \X{x 1 -x f, (3.6) 

with xq an arbitrary offset and A, ui, 2tt an incommensu- 
rate triplet. (We take A = 2, uj = 2n/\/5.) From studies 
of the d-dimensional quantum kicked rotator it is known 
that such a simple potential, which is linear in d — 1 di- 
mensions and nonlinear in one single dimension, provides 
sufficient randomness for localization^ 



IV. LOCALIZATION IN THE QUANTUM HALL 
EFFECT 



A. Numerical simulation 

We base our numerical simulation on the ID strobo- 
scopic model with two incommensurate frequencies of 
Sec. IIIII We introduce a Bloch number q and seek the 
time dependence of the state ip(pi,t) = e~ lqpi x q {pi, t). 
The state Xq(pii~t) 1& a 27r-periodic function of p\, so it is 
a superposition of a discrete set of eigenstates e~ lmpi of 
x\. For numerical purposes this infinite set is truncated 
to M states, m £ {1, 2, . . . M}, with periodic boundary 
conditions at the end points. 

Fourier transformation from eigenstates of x\, with 
eigenvalue m, to eigenstates of pi, with eigenvalue 
2im/M, amounts to multiplication with the unitary ma- 
trix 



M -l/2 e 2rinm/M i 71, 771 G {1, 2, . . . M}. (4.1) 



Calculation of the state Xq{ x ii t), f° r t an integer multiple 
of t = 1, requires 2t Fourier transformations, 



't-i 



Xq(xi,t)= Y[F q (t')]xq(xi,0), 



(4.2) 



U'=0 
M 



[FMnm = J2 U k^~ iH0{2Wk/MM+a) Uk m e- iVl(m+q) . 

(4.3) 



fe=i 



These operations can be carried out with high efficiency 
using the fast-Fourier-transform algorithm^ 
As initial state we choose 



Xq(xi,0) = 8 X0 . X1 



-i<j>o/2 



cos(0o/2) [ I 



^/ 2 sin(0 o /2) (I 



(4.4) 
(4.5) 



spatially localized at x\ = xo = M/2 (for even M). The 
angles 4>o > of the initial spin direction are chosen ran- 
domly on the unit sphere. 



B. Localization-delocalization transition 

To search for localization we calculate the expectation 
value 

M 

(( Xl {t)-x ) 2 ) = J2(™-M/2) 2 \ X q(xi=m,t)\ 2 , (4.6) 



and obtain the mean squared displacement 



A 2 (t) = ((x 1 (t)-x Q y) 



(4.7) 
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FIG. 2: Time-dependent diffusion coefficient (£g| at t = 3000 
as a function of p for M = 1024. The different values of K 
range from weak disorder (K — 4) to strong disorder (K — 
0.8). The peaks signal a localization-delocalization transition. 
Compared to the three quantum Hall transitions in a clean 
system (indicated by arrows), the two outer transitions are 
displaced inwards by disorder, while the central transition is 
split. The splitting of the two central peaks becomes larger 
and larger with increasing disorder, until they merge with the 
outer peaks. 



by averaging over some 10 2 -10 3 values of the random 
parameters a,q £ {0,1}. We fix K = 2, f3 = 0.8 and 
vary the parameter p. 

The system is localized if the time-dependent diffusion 
coefficient 



D(t) 



A 2 (i) 



(4.8) 



vanishes in the large- t limit. Delocalization with diffu- 
sive propagation corresponds to a non-zero large-time 
limit of D(t). The quantum Hall phase transition is a 
localization-delocalization transition, so we would expect 
a peak in D(t) as a function of p at the critical points p c 
where the topological invariant switches from one value 
to another. In a clean system these values are p c = 0, ±2, 
see Eq. (j2~2|) . 

The data in Fig. [2] shows that disorder has two effects: 
It shifts the outer transitions inwards and splits the cen- 
tral transition, resulting in a total of four peaks. We will 
demonstrate in Sec. [V] that these are topological phase 
transitions, by calculating the topological invariant - 
which as we can see in Fig. [3] switches at each of the 
transitions. 



C. Scaling and critical exponent 

In the single-parameter scaling theory of localization 
all microscopic parameters enter only through a single 
length scale £ (the localization length) and the associ- 
ated energy scale 6% = (£ d p c ) 1 (being the mean level 
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FIG. 3: Left axis: Time-dependent diffusion coefficient for 
K — 2 (solid curve, same data as in Fig. [2]), showing the 
four localization-delocalization transitions. Right axis: Four- 
terminal Hall conductance Gh (data points) and topological 
invariant I (dashed curve), calculated in Sec. [Vj to demon- 
strate that these are topological phase transitions. 



spacing in a c?-dimcnsional box of size £, obtained from 
the density of states p c at the critical energy) The 
corresponding scaling law for dynamical localization has 



D(t) = ^- d F(C d t) 



(4.9) 



in the large-time limit near the critical point p c . The 
localization length £ diverges as a power law with critical 
exponent v on approaching the transition, 



£ oc p c \~ 



(4.10) 



The limiting behavior of the function F(z) is F(z) oc 
1/z for z — > oo and F(z) oc z 2 ^- 1 for z -> 0. The 
first limit ensures that the mean squared displacement 
A 2 = tD(t) — > £ 2 becomes time independent in the limit 
t — > oo at fixed fj, — fi c . The second limit ensures that, if 
we send /j, —> fx c at fixed t, the diffusion coefficient D(t) — > 
f.2/d-i t cnc ! s to a finite value. For d = 2, this value 
is also time independent, which implies regular diffusion 
(D = constant) at criticality in two dimensions. 

We have performed a finite-time scaling analysis of 
D(t), similar to Refs. HU to obtain the localization 
length £ and extract the value of the critical exponent. 
(See App. [B] for details.) We considered times up to 
t = 1.3 ■ 10 6 for system size M = 2 13 = 8192. In Fig. H 
we show both the unsealed and the scaled data. For the 
two independent transitions we find v = 2.576 ± 0.03 at 
p c = 0.387 and v = 2.565 ± 0.03 at p c = 1.903. Both 
results agree with i^qhe = 2.593, the critical exponent 
for the quantum Hall phase transition^ 



V. HALL CONDUCTANCE AND 
TOPOLOGICAL INVARIANT 

The anomalous quantum Hall effect in the absence 
of disorder (V = 0) is characterized by the topological 
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FIG. 4: Top panel: Time-dependent diffusion coefficient for 
different times as a function of /i. The curves are a least- 
squares fit, used to extract the localization length and 
the critical exponent (see App. [Bj. In the lower panel time 
is rescaled, to test the scaling form D(t,n) = F(t/^ 2 ) [see 
Eq. (|4.9[) ]. The open data points do not fully collapse onto a 
single scaling curve, due to finite-time corrections to scaling. 
The filled data points include the leading-order correction (see 
App. El). 

invariant^ 3 - 

i r , r , \du(p) du(p)] , 

X=-—\ d Pl d P2 — — x — W -u(p), 
47T J -it J-* [ dpi dp 2 

(5.1) 

with u = u/\u\. This socallcd Skyrmion number does 
not apply for nonzero disorder potential, when momen- 
tum p is no longer a good quantum number. 

We calculate the topological invariant for nonzero V 
from the winding number of the reflection matrix r(</>) in 
a cylinder geometry 38 ! 39 

X =--!/" d0-^lnDetr(0), (5.2) 
2m J Q dp 

where $ = 4>h/e is the flux enclosed by the cylinder and 
r((j>) is evaluated at e = 0. (We explain in App. [C] how 
to construct the quasi-energy dependent reflection ma- 
trix from the Floquet operator . 40 ' 41 ) Since this is a 2D 
system, the sizes Mx M for which we can calculate I are 
much smaller than in the ID reduction used to calculate 
D(t). 

The results in Fig. [5] are for M x M = 40 x 40. This 
is data for a single sample (q = 0), at fixed K = 2 as a 
function of /?, fx. The disorder-averaged /^-dependence of 
X is plotted in Fig. [3] (dashed curve, for (3 = 0.8). 

Comparing with the phase boundaries (|2.2[) for the 
clean system (V = 0), we see that disorder introduces 



topologically trivial regions along clean phase bound- 
aries. In the disordered system transitions between two 
different topologically nontrivial phases (with I = ±1) 
go via a topologically trivial region (1 = 0). A similarly 
disruptive effect of disorder (but with a metallic gapless 
region replacing the topologically trivial phase) has been 
observed in computer simulations of the quantum spin 
Hall effect^ 

We have also calculated the Hall conductance Gh, 
which unlike the topological invariant is a directly mea- 
surable quantity. The results shown in Fig. [3] (data 
points) were obtained in a single four-terminal sample 
of dimensions M x M = 70 x 70, directly from the scat- 
tering matrix expression for the Hall conductance^ 3 - The 
Hall plateaus are at the values expected from the topo- 
logical invariant, Gh ~Ix e 2 /h, with deviations from 
exact quantization due to the relatively small size of the 
2D system. 

VI. DISCUSSION 

We have shown how the quantum Hall effect can be 
modeled in a ID dynamical system, by using a pair of 
incommensurate driving frequencies to simulate the ef- 
fect of a second spatial dimension. This ID stroboscopic 
model could become a competitive alternative to the 
2D network model for numerical studies of the quantum 
Hall phase transition^ similarly to how the ID quantum 
kicked rotator is an alternative to the 3D Anderson model 
of the metal-insulator transition^ 

Since quantum kicked rotators can be realized using 
cold atoms ; 10 ' 11 ' 37 ' 44 the stroboscopic model might also 
provide a way to study the quantum Hall effect us- 
ing atomic matter waves. Cold atoms represent clean 
and controllable experimental quantum systems, owing 
to the ability to tune interaction strengths and external 
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potentials^ Due to the absence of impurities they have 
long phase coherence times, so their quantum dynamics 
can be followed over long time scales. These properties 
make cold atoms ideally suited for the experimental study 
of quantum phase transitions. 

There is a particular need for a new physical system 
to investigate the quantum Hall phase transition, because 
currently the theory disagrees with semiconductor exper- 
iments on the value of the critical exponent.— This might 
be an effect of Coulomb interactions between the elec- 
trons in a semiconductor, and a system with controllable 
interactions could shed light on this question. 

For cold atomic gases prepared in a magneto-optical 
trap a quasi-periodically modulated ID standing wave, 
created by two overlapping laser beams, simulates the 
quasi-periodic driving of the kicked rotator4i The mo- 
mentum distribution is accessible through an absorption 
measurement, following the release of the atomic gas from 
the trap4^ Since in the kicked rotator momentum plays 
the role of coordinate, in this way the diffusion coef- 
ficient could be measured and the critical exponent of 
the metal-insulator transition was obtained from its time 
dependence 1 10 i 11 i 37 

To realize the stroboscopic model of the quantum Hall 
effect, a controllable spin-f jl degree of freedom is needed. 
Hyperfine levels in alkali or earth alkali atoms can be used 
for that purpose,— and arbitrary rotations of this pseu- 
dospin have been demonstrated in Cs4^ Two overlapping 
standing waves would produce a purely sinusoidal kicking 
potential (corresponding to T(u) = 1), while for flat spin 
bands higher harmonics are desirable. Fortunately, the 
topological nature of the phase transition ensures that 
there is considerable freedom in the choice of the poten- 
tials. 

Continuing on the path of dimensional reduction pro- 
posed here, it is conceivable that the hypothetical 4D 
quantum Hall effect 47 might also be realized in the labo- 
ratory, by adding two more incommensurate driving fre- 
quencies. 
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Appendix A: Tight-binding representation 

To gain further insight into the stroboscopic model, we 
give a tight-binding representation. This will motivate 
the specific form (12.51) for the function T(w) and it will 
also guide us in the choice (|2.6[) for the scalar potential 



V(x). The derivation follows the same steps as for the 
quantum kicked rotator^^ 

Including the spin degree of freedom (s = ±1), we 
denote the coordinate basis states by \m,s), such that 
Xi\m, s) = rrii\m,s) and cr z \m, s) = s\m, s). The two 
states \a±) defined by 



T q \a+) = e I£ |a+), 



\a^)=e iH °\a 4 



e ie - lV "\a+) 



(Ala) 
(Alb) 



are evaluated just after and just before the kick. [We have 
abbreviated V q = V(id p + q).] Both states evolve with a 
phase factor e~ le in one period r = 1. The tight-binding 
representation is expressed in terms of the average 



|6) = i(|a + ) + |a_)). 
The Hermitian operator 



(A2) 



f — e iHa 1 
W = i -rr = - taxi\\uT{u)\u ■ er (A3) 



allows to relate \b) to \a±) separately, 
1 , I 



\b) 



(A4) 



1 + iW 1 ' 1-iW 
Substitution into Eq. (|Alb[) gives 

(I + iW)\b) = e ie - iV «(l-iW)\b) (A5a) 

(A5b) 



1 _ pie-iVq 

i - — \b) = W\b) 



1 



oie—iV a 



t^[{e-V q )/2}\b)=W\b). 



(A5c) 



In coordinate representation this gives the tight- 
binding equations 



W{m+q)-\s 



¥ m = 0, (A6) 



with hopping matrix elements 

W„' = (m,s\W(p)\m + n,s'), (A7a) 
W(p) = -tim\±uT(u)]u- a. (A7b) 

The tangent term in Eq. (|A6|) provides a pseudo- 
random on-site potential, provided that V(m + q) 
changes from site to site in a way which is incommen- 
surate with the periodicity it of the tangent. This is why 
a simple polynomial V(x) suffices to produce the local- 
izing effect of a disorder potential^ 

The role of the Bloch vector q is to provide different 
realizations of the disorder potential, so that a disorder 
average is effectively an average of q over the Brillouin 
zone. The strength of the disorder potential is varied by 
varying the parameter K, which determines the relative 
magnitude of kinetic and potential energies: small K 
corresponds to strong disorder. 
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From Eq. (|A7[) we see that different choices for T(u) 
lead to different hopping matrix elements, leaving the 
on-site disorder unaffected. The arctangent form in Eq. 
(|2.5[) has the simplifying effect of excluding hopping be- 
tween sites that are not nearest neighbors. For this choice 
— tan[iuT(u)] = 1 the hopping matrix elements are given 

by 

W n = 2irKl3Li<T z 5 nu0 Sn 2 ,o 

+ nK(±ia y - p(T z )5 ni fiS n2t ± 1 

+ nK(±ia x - l3a z )S„ u ±iS n2fi . (A8) 



Appendix B: Finite-time scaling 

Following Refs. @H3, we extract the critical exponent 
v from finite-time numerical data by fitting the diffusion 
coefficient (or, more conveniently, its logarithm) to the 
scaling law D(t) = F(t/^ 2 ). For finite t the diffusion co- 
efficient is an analytic function of fi. In view of Eq. (|4.f 0[) 
the variable {t/^ 2 ) 1 / 2lJ = t x l* v u is an analytic function of 
/i, vanishing at fi c . 

We therefore have the two power series 

JVi k 
lnD(t) =lnD c + ^4 1) (t 1 / 2v u) + c n r\ (BI) 

k=l 

N 2 

U = /i-/i c + ^ci 2) ( A1 - Alc ) fc . (B2) 

k=2 

The term c$t~ v , with y > 0, accounts for finite-time 
corrections to single-parameter scaling at the transition 
point. We put c\ = 0, c-i < 0, to ensure that D(t) as 
a function of \x has a maximum at /i c . Wc then choose 
integers N\,Ni and fit the parameters D Cl v,co,y with 
c fc (' G {1, 2}, 2 < k < Ni) to the t and /i-dependence 
of D(t), for a given ID system size M. 

For the transition around /i = 0.38 wc took times 
t = 1.2 • 10 4 , 3.3 • 10 4 , 8.3 • 10 4 , 2.1 • 10 5 , 5.2 • 10 5 , and 
1.3 • 10 6 , with M = 2 13 = 8192. We averaged over 1000 
samples. The quality of the fit is quantified by the chi- 
square- value per degree of freedom (x 2 /ndf). We system- 
atically increased N\,N2 until we arrived at x 2 / n df ~ 1. 
Only the leading order term in Eq. (|B2[) was needed for 
a good fit, so we simply took u = fj, — /i c . The expan- 
sion (|Bip did need higher order terms, up to Ni = 6. 
We thus obtained v = 2.576 ± 0.03 at /i c = 0.387 with 
X 2 /ndf = 1.2. A similar analysis was performed for the 
outer peak in fig. O resulting in v = 2.565 ± 0.03 at 
At c = 1-903 with x 2 /ndf = 1.01. 



Appendix C: Scattering matrix from Floquet 
operator 

To calculate the topological invariant (15.21) wc need the 
reflection matrix r((j)) in a cylinder geometry, at quasi- 



1 ► mi 

FIG. 6: Truncation of the lattice used to construct a four- 
terminal scattering matrix, as described in the text. 



energy e = as a function of the flux $ = (j>h/ e en- 
closed by the cylinder. This can be obtained from a four- 
terminal scattering matrix 5, which relates the wave am- 
plitudes of incoming and outgoing states at the four edges 
of an M x M square lattice of sites (3:1,22) = (^1,7722), 
mi = 1,2,. . . M. The dimensionality of S is 8M x 8M, 
with the factor of 8 accounting for four terminals and a 
twofold spin degree of freedom. The Floquet operator T q 
is a 2M 2 x 2M 2 matrix describing the stroboscopic time 
evolution of states on the 2D lattice. (We do not make 
the dimensional reduction to ID for this calculation.) 

When the square is folded into a cylinder, incoming 
and outgoing states at the left and right edge are re- 
lated by a phase factor e 1 ^. This relation can be used 
to reduce the four-terminal scattering matrix to a two- 
terminal scattering matrix S(4>) (which now has dimen- 
sion AM x AM). The reflection matrix r{4>) is a 1M x 2M 
subblock of S((f>), relating incoming and outgoing states 
at the lower edge. We refer to Ref.[39|for a computation- 
ally efficient way to carry out this general procedure. 

What we discuss in this Appendix is how to obtain S 
from T q . We are faced with the complication that the 
truncation of the coordinates to a finite range M intro- 
duces spurious hopping matrix elements that couple sites 
near opposite edges (typically within 5-10 sites from the 
edge). We cannot directly delete these matrix elements 
from the Floquet matrix without losing unitarity. 

Our solution (illustrated in Fig. [5]) is to start from a 
larger M' x M' system (red square), with Floquet ma- 
trix F' q = e- lH '»e- lV «. We then go back to the M x M 
system (green square) by deleting rows and columns in 
the coordinate representation of H' M- Ho and V' 1— ^ V q . 
The resulting Floquet matrix T q = e~ lH ° e^ 1 ^ remains 
unitary. By choosing M' sufficiently larger than M (typ- 
ically M' = M + 10 suffices), we effectively eliminate the 
spurious hopping matrix elements. 

For a four-terminal scattering matrix we introduce ab- 
sorbing terminals at the four edges of the M x M lattice. 
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The 8M x 2M 2 matrix P projects onto these terminals, 

C 1 if mi G {1,M}, 

-Pmm' = <5 S a'* mi m'/m a mi X <M if TO 2 € {1,M}, (CI) 

I otherwise. 

The e-dependent scattering matrix S is obtained from 
the Floquet matrix T q through the formul a 40 ' 41 

S = P [l - e le Jq(l - P T P)] _1 e ie F q P T , (C2) 

where the superscript T indicates the transpose of the 
matrix. The quasi-energy s is set to zero for the calcula- 
tion of the topological invariant (j5.1[) . The integral over 
</> is evaluated analytically^ by contour integration over 



complex z = e 1 ^ . Results for M — 40 are shown in Figs. 
E] and EJ 

To calculate the Hall conductance Gh we directly use 
the four-terminal scattering matrix S, without rolling up 
the system into a cylinder. The geometry is still that of 
Fig. |6l but the four absorbing terminals are point con- 
tacts, covering a single site at the center of each edge. 
The dimensionality of the scattering matrix, including 
spin, is thus 8x8. A current /13 flows from terminal 1 
to terminal 3 and the voltage V24 is measured between 
terminals 2 and 4 (which draw no current). The Hall 
conductance Gh = /13/V24 is obtained from the scatter- 
ing matrix elements using Biittikcr's formulas^ Results 
for M = 70 are shown in Fig. [3] 
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